Bull. Spec. Astrophys. Obs., 2011, 66, 199-206 



© Special Astrophysical Observatory of the Russian AS, 2011 



Statistics of WMAP ILC Map Temperature Fluctuations 
Towards Distant Radio Galaxies 
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For 2442 galaxies of the catalog, compiled based on the NED, SDSS, and CATS survey data 
with redshifts z > 0.3 we conducted an analysis of the amplitude of temperature fluctuations in 
the cosmic microwave background (CMB) in the points, corresponding to the direction to these 
objects. To this end, we used the ILC map from the WMAP mission seven-year data release. We 
have estimated the dipole component of the background and tested the hypothesis of Kashlinsky 
on the existence of a "dark bulk flow" , determined for the estimated dipole component of the CMB 
WMAP by the value of the CMB anisotropy in the direction to the clusters of galaxies. We show 
that the amplitude of this dipole T max — 0.012 mK is located within the a interval, estimated by 
Monte Carlo simulations for the Gaussian fluctuations of the CMB signal in the ACDM model. 
The low amplitude of the dipole indicates that it is impossible to confirm this hypothesis from the 
WMAP data. In addition, we studied the statistics of the fluctuation amplitude of the microwave 
signal in the direction to radio galaxies. A weakening of the absolute value of the mean signal in 
the corresponding fields was discovered. 
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1. INTRODUCTION 

The studies of properties of the large-scale structure, 
reflected in the statistics of the Cosmic Microwave 
Background (CMB), are of particular interest in con- 
nection with the advent of new high-accuracy obser- 
vations in the millimeter and submillimctcr ranges 
[3 [21 [3J H]. Using the data on galaxy clusters and 
the magnitude of the CMB signal fluctuations in the 
direction to their centers, where the signal fluctu- 
ations are caused by the Sunyaev-Zel'dovich effect 
(SZ) [5], we can estimate the cosmological parame- 
ters (the Hubble constant Hq, and the density pa- 
rameter fio) 13 ■ I R addition, one of the components 
of the SZ effect — the kinematic SZ effect, caused by 
the proper motion of the cluster as a whole with re- 
spect to the CMB reference system — opens the pos- 
sibility of studying the large-scale fluxes of matter. 
In particular, [8] describe the study of large-scale 
"dark bulk flows" in motion of clusters of galaxies, 
based on the WMAP5 data and the catalog, com- 
piled from the X-ray data and containing more than 
1000 clusters. The authors have calculated the dipole 
from the CMB pixel values. The dipole is determined, 
in their opinion, by the kinematic SZ effect; it reveals 
a constant velocity, at least on the scale of 800 Mpc. 



Note that despite the low WMAP sensitivity to 
the SZ effect due to the relatively poor resolution of 
the final map of the isolated CMB signal (the lim- 
iting spherical harmonics — a multipole — in the first 
WMAP releases for the ILC CMB map there is 
imax ~ 100, or, in terms of resolution, ~55' and in 
the seven-year release it is l max ~ 150 (8 ~35'), while 
the SZ effect for most of galaxy clusters is strongly 
manifested on the small scale: 9 <10'), the discussed 
effect is of great interest in terms of its manifestation 
in the ongoing Planck missio 

Also note that despite the limitations on reso- 
lution (channel W has a better resolution of 12.6' 
with the frequency of 94 GHz) and on the sensitiv- 
ity of WMAP maps, various attempts to measure 
the SZ effect statistically have been made both in 
the population averaging ("stacking") (see, e.g., re- 
cent papers [H[T0]), and directly for individual bright 
clusters, such as Coma pi)], where the effect was 
observed. Diego and Partridge [9] used three chan- 
nels: Q (43 GHz, FWHM = 0.51°), V (61 GHz, 
0.35°) and W (94 GHz, 0.22°) and obtained aver- 
aged profiles in the direction to known clusters. In 



1 http://www.rssd.esa.int/Planck/ 
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particular, besides the presence of the desired signal 
and the growth of its amplitude with frequency, the 
authors found that the signal's value is lower than ex- 
pected from the X-ray data, the fact that they linked 
with the presence of a point source in the cluster. The 
WMAP team [10] has conducted a similar study, us- 
ing two channels: V and W. They also applied the av- 
eraging of different regions in the direction to known 
clusters with X-ray emission and evidenced the ef- 
fect. In addition, they discovered that the effect can 
be traced down to the scales of 6 — 1.05°. The sig- 
nal on the scale of 8 = 0.58° is considered as true, 
and that on the scale of = 1.05° — as a statistical 
fluctuation. At that, the authors [10] have concluded 
that the apparent effect for the averaged sources is 
consistent with that, expected from the X-ray obser- 
vations. Let us emphasize that the evaluations are 
made for the source, producing the SZ effect, and av- 
eraged over the sky. Recent experiments, such as the 
the one conducted on the Atacama Cosmology Tele- 
scope (ACT) [IT] and Planck mission [12] already al- 
low to isolate and see this effect directly. 

Despite the fact that for most clusters the SZ ef- 
fect is not directly visible in the WMAP data, it is 
estimated statistically. We believe that the approach 
proposed in [7J [S] deserves special attention. We de- 
velop it in this paper, using the amplitude of temper- 
ature fluctuations in the points, corresponding to the 
direction to distant radio galaxies from our survey for 
the analysis. To detect the possible "dark bulk flow" 
we need to evaluate the signal towards distant radio 
galaxies. A useful moment is the fact that the CMB 
signal in the ILC map of fluctuations is isolated using 
the channels, in which the SZ effect, if existing, will 
have a negative signal. Fortunately, the angular scale 
of the SZ effect, discovered by the WMAP team [TU] 
corresponds to the resolution of the ILC map of seven- 
year observations. 

In the standard scheme of galaxy formation, the 
radio source lights up as a result of a merger of galax- 
ies, and formation of an accretion disk and jets, ob- 
served at radio and other wavelengths. As a rule (see 
reviews in [13[ 114]). the most powerful radio galax- 
ies, visible at large redshifts are identified with gi- 
ant elliptical galaxies, which are mainly the central 
galaxies of clusters and are formed by merging. This 
property (radio emission) can be used to find dis- 
tant clusters and protoclusters of galaxies. For ex- 
ample, [15] presents the results of a program, con- 
ducted with the ESO VLT telescope, searching for 
the emerging clusters of galaxies near the powerful 
radio galaxies at redshifts 2 < z < 5.2 with the ra- 
dio fluxes i?2.7GHz > 10 33 erg s _1 Hz -1 sr _1 . The 
authors have selected 150 objects and examined the 
fields around nine of them. In the fields the authors 



have selected the galaxies, emitting in Lya (the so- 
called Lya-emitters) , the redshifts of which were mea- 
sured. Using the data on spatial density of objects it 
was concluded whether they belong to protoclusters. 
The size of such protoclusters is at least 1.75 Mpc. It 
was shown that 75% of radio galaxies with z > 2 are 
associated with protoclusters. Hence, we can estimate 
that approximately 3xl0~ 8 of the forming clusters 
lay in the interval of 2 < z < 5.2 in the comoving 
unit volume with a side of 1 Mpc with an active radio 
source. But it is likely that within a given range of 
redshifts the number of protoclusters may be higher, 
since an active radio source may simply not be ob- 
served. 

Using the survey of remote (for example, z > 0.3) 
radio galaxies, we can track the locations of clus- 
ters and protoclusters of galaxies, which in turn 
may affect the fluctuations of the microwave back- 
ground. A sample of radio galaxies, prepared by us 
earlier [THl H3 HB] is very suitable to this end, and 
we use it in our study. 

2. DATA ANALYSIS 

2.1. Input Arrays 

A sample of radio galaxies with 

z > ' 0.3 i EH i i 1 

was constructed using the tools and archives of 
the following databases: NEE@, CATS^I [21] [22], 
SDSS0 23 with the aim to use it in various statisti- 
cal and cosmological tests (T3l [2H [25] , which require 
an analysis of a large sample of objects of the same 
nature. The NED database was used to construct 
the primary list of objects. We selected from it the 
objects with required parameters, mainly based on 
redshift (z > 0.3) and morphological properties of 
radio galaxies. The initial catalog contained 3364 
objects. Such a sample of galaxies is contaminated 
with objects with incomplete information, or objects 
with other properties. Therefore, special attention 
was paid to clean the original sample from the odd 
sources: galaxies (1) with redshifts, determined by 
the photometric method; (2) those with quasar 
properties from the available literature data. The 
final catalog contains 2442 sources with spectroscopic 
redshifts, the photometric values and flux densities 
in the radio range, the sizes of radio sources, as well 
as radio spectral indices, which were calculated based 
on the results of cross-identification with the radio 
catalogs, stored in the CATS, in the frequency range 



2 http://nedwww.ipac.caltech.edu 

3 http://cats.sao.ru 

4 http://www.sdss.org 
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from 325 MHz to 30 GHz. The positions of 2442 
galaxies on the sphere are demonstrated in Fig. 1. 

As the CMB map, we used the map of 
the distribution of the cosmic microwave back- 
ground anisotropy, reconstructed based on the multi- 
frequency observations via the technique of internal 
linear combination (ILC) of the background compo- 
nents p] in the observational program of the WMAP 
space mission^ (Wilkinson Microwave Anisotropy 
Probe). To construct the WMAP map, the data 
in five bands were used: 23 GHz (K-band), 33 GHz 
(Ka-band), 41 GHz (Q-band), 61 GHz (V-band), and 
94 GHz (W-band). The ILC map contains informa- 
tion on the distribution of the microwave background 
for not very high harmonics (£ < 150). 

Nevertheless, it is currently the most tested and 
used CMB map. For the data analysis, we use the dis- 
tribution of the signal AT (#,</>), describing the tem- 
perature anisotropy on the sphere, with two limit- 
ing values (by angular resolution) of the multipoles 
e max < 150 (9 >36') and £ max < 20 (9 >260') 
according to the decomposition into spherical 
harmonics: 



oo m—l 



AT(0, 4>) = atmYt 



Cm \ 



i=2 m=-l 

where the spherical harmonics is 



Yt 



tin \ 



(21+1) (l-m)! 
4-7T (£ + m)\ 



P £ m (x)e 



x = cos 9 Pp(x) are the associated Legendre poly- 
nomials, £ and m are the multipole number and its 
modes, respectively. For the expansion of spherical 
harmonics, we used the GLESP packaged [26 . 

2.2. Signal Histograms in CMB Pixels 

To analyze the Gaussianity of the distribution in the 
first order and compare it with the model, we took 
the amplitude of CMB signal fluctuations in the map 
pixels, onto which the radio galaxies are projected 
(Fig. 2), and constructed the histograms of its distri- 
bution (Fig. 3). One of the rules of galaxy selection 
in our sample was the presence of a spectroscopic 
redshift, measured, as a rule, in objects outside the 
Galactic plane, hence the lack of sources in the cen- 
tral band of Fig. 2 (the Galactic plane) is a selection 
effect. The pixel size of maps, used for the analysis, 
at resolutions of £ max < 150 and £ m ax < 20 amounts 
to 36'x36' and 260' x 260', respectively. The analysis 



5 http://lambda.gsfc.nasa.gov 

6 http://www.glesp.nbi.dk 



was carried out for maps with £ max < 150, since the 
ILC data in the WMAP7 release are presented with 
a better resolution. 

Figure 3 shows the corresponding histograms for 
the CMB fluctuation amplitude distributions of the 
ILC map in pixels, corresponding to the direction to 
the catalog radio galaxies. The histograms present a 
lcr-spread, calculated by two methods on the pixel 
statistics: (1) for the given 100 realizations of the mi- 
crowave background maps in the ACDM cosmologi- 
cal model with homogeneous and isotropic Gaussian 
random fields, resulting in the corresponding CMB 
fluctuations, and (2) for 100 realizations of a random 
arrangement of 2442 points in the ILC WMAP7 map 
itself. The value of the histogram bean is 0.02 mK. 

Note the particularities of the distributions shown 
in Fig. 3. All the signal variation distributions in the 
studied pixels are close to normal. The distribution 
maximum N is located strictly in the region of zero 
amplitude fluctuations. Its value exceeds the mean ex- 
pected level by more than 3er for the maps with a res- 



olution of . 



< 150 and about 2a with L 



< 20 



for the spread estimates, calculated using the Monte 
Carlo method for Gaussian background fluctuations 
in the ACDM cosmology (left graphs of the figure), 
and over 2a for the estimates, obtained in randomly 
selected pixels of the ILC map. The estimates are pre- 
sented in the Table. 

In addition, for both resolutions £ = 150, and 
£ 20, the ratio (T maa: y/T m i n ) > (T ma x/T m i n y , 
indicating a lower dispersion of signal variations in 
the map than that, expected in the ACDM model. 

2.3. Dipole Estimation 

Using the value of temperature fluctuations in pixels 
in the map with a resolution of I = 150, applying the 
least-squares method we estimated the corresponding 
dipole in the form 

T(l, b) = T x cos(Z) cos(b) + T y sin(Z) cos(fe) + T z sin(fo), 

where (I, b) are the galactic coordinates: longitude 
and latitude, T(Z, b) is the value of CMB temperature 
fluctuations in mJy, taken in the pixels corresponding 
to the positions of radio galaxies, and (T x , T y ,T z ) are 
the dipole components, which are found to be equal 
to (0.0116, 0.0036, 0.0026), respectively. 

Hence, we obtain the following for the extrema: 



T, 



±JT*- 
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b = ± arctan^/JTj + T y 2 ), 
I = arctan(T y /Tc), 
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Figure 1: Position of selected radio sources on the celestial sphere in galactic coordinates. White circles mark 
the SDSS objects, gray crosses — the remaining sources. 




Figure 2: Position of selected galaxies (black circles) in the CMB WMAP map with a resolution of i max = 20 
in galactic coordinates. Dark spots in the map correspond to the cold signal, while light spots — to the hot signal. 



or in the maximum T max = 0.0124 mK with co- 
ordinates (lmax,b ma x) = (17° .307, 1 1° .844) , and in 
the minimum T TO j„ = — 0.0124mK with coordinates 
{lmin,b mi n) = (197°.307, -11°.844). The dipole's po- 
sition on the sphere is shown in Fig. 4. 

Despite the fact that the estimated dipole is vi- 
sually "drawn" towards the Galactic center, the sig- 
nificance of its position, as well as the reality of its 
existence seems doubtful in view of its small ampli- 
tude: T max = 0.012 mK. This value lays within the o- 
scatter from the mean, i.e. within the noise, estimated 
by modelling of 50 CMB signal realizations in the 
ACDM model, and determining the parameters of the 



dipole in each of the realizations for the coordinates of 
the radio galaxy catalog. The modelling results yield 
an average estimate of the parameters (T x ,T y ,T z ) = 
(0.003 ± 0.016, 0.001 ± 0.018, -0.002 ± 0.010). 

3. RESULTS 

We investigated the properties of the CMB 
WMAP7 signal in the fields of distant (z > 0.3) radio 
galaxies. In general, the signal distribution over the 
areas corresponds to normal. Moreover, in this dis- 
tribution the histogram amplitude is higher than ex- 
pected (more than by a cr-scatter) both for the ACDM 
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Figure 3: The histograms of the distribution of the microwave background signal value (the thick solid line) 
in the ILC map, measured in pixels, corresponding to the direction to radio galaxies. Above: histograms for 
the pixel map resolution of l max < 150. Bottom: histograms for the pixel map resolution of £ max < 20. Left: 
the dotted and dashed lines show the levels of la and 2a-spreads, respectively, calculated from the data of 100 
realizations of the ACDM cosmological model. Right: The dotted and dashed lines show the levels of la and 
2(7 -spreads, respectively, calculated from the data of 100 realizations of a random distribution of points in the 
ILC WMAP seven-year map. 




Figure 4: The position of the dipole on the sphere, determined with the least-squares method from the CMB 
pixel values in the locations of distant radio galaxies, in galactic coordinates. 
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Table 1: Table. Parameters of normal distribu- 
tion (the amplitude N and the width parameter 
sq = 0o.5(81n2)~ 1 / 2 J, corresponding to the distribu- 
tion of the CMB temperature fluctuation amplitude 
in the pixels of galaxy positions for maps with res- 
olutions I = 150 and I = 20. The estimates of the 
distribution of temperature minima and maxima at 
the a-scatter are obtained when modeling the back- 
ground fluctuations in 100 realizations of a random 
signal in the KCDM cosmological model (marked as 
ledm), and in the simulations of 100 realizations of 
a random scatter of galaxy locations in the ILC map 
(noted as ilc). 



Value 


Amplitude 


so, 


Distribution 


N 


mK 


T(£ = 150) 


319.2 


0.060 


T min {l= 150) icdm 


232.7 


0.072 


T max (l= 150) lcdm 


271.0 


0.082 


T m in{@ — 150) 1 C 


276.1 


0.064 


T max (t = 150) 4 c 


304.4 


0.068 


T(£ = 20) 


552.6 


0.035 


T mm (£ = 2Q) lcdm 


362.4 


0.041 


T max (£ = 20) lcdm 


501.1 


0.049 


Tmini,^ — 20) 


438.8 


0.041 


Tmax [£ — 20) 


474.5 


0.043 



cosmological model and for the random positioning 
of objects on the ILC WMAP map. That is, the pix- 
els, corresponding to the direction to distant radio 
galaxies, are dominated by small (approaching zero) 
fluctuation values, which leads to an increase in the 
histogram amplitude above the expected value. Ab- 
sence of a shift in the position of the histogram peak 
indicates either the absence of any signal in these pix- 
els, or a compensation of the SZ effect by the radio 
emission of the galaxy (see, e.g., [9]). Theoretically, 
one could expect a shift either towards the negative 
signal in the case of the SZ effect, or towards the 
positive signal in the case of existence of a residuary 
signal from the radio source after the component sep- 
aration. In the presence of both factors we would see 
an increase in the width of the histogram, compared 
with the simple models of the background perturba- 
tion. However, an inverse effect is observed: the re- 
sulting CMB WMAP7 signal distribution, measured 
in the regions of radio galaxy positions on the sphere 
is in fact narrower than the one, expected for random 
Gaussian fluctuations, and than the one, observed in 
the ILC map on the average. In addition, the variance 
of the so distribution is smaller than expected in mod- 
els, which also indicates an increased number of pix- 
els with zero signal in the fields of radio galaxies. And 
if the correlation {T max /T mm ) lcdm > {T max /T mm ) Uc 



can be explained by an underestimated value of the 
amplitude of the ILC quadrupole, then the cause of 
the effect of signal attenuation in the fields of radio 
galaxies is as yet unclear. 

We as well tried to test the effect of existence of 
a dipole in the CMB data in the regions of galaxy 
clusters, discovered in [3 [8] . It is determined, as the 
authors suppose, by the kinematic SZ effect, and asso- 
ciated with the "dark bulk flow" of matter. To verify 
this phenomenon, we constructed a dipole using the 
values in the pixels of galaxy positions, and found 
that this dipole's amplitude is below the a noise level 
of model variations. 

Note that the size of the ILC map pixel in general 
does not allow to valuably investigate the SZ effect. In 
addition, it is clear that the nonzero dipole can almost 
always be defined for a finite sample of pixels with 
nonzero values, what the models are demonstrating. 
Although we do not rule out the possibility that for 
the nearby clusters of galaxies the background fluctu- 
ations can show the existence of a common movement, 
distinct from the known CMB kinematic dipole, but 
for distant objects, as in our case, the Hubble flow will 
be dominating. We shall be able to check this effect 
after the publication of maps of the Planck mission. 
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